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Abstract: It is pointed out that there exist a few problems to be overcome toward an 
observable sub-eV QCD axion in superstring compactification. We give a general expression 
for the axion decay constant. For a large domain wall number N^Wi the axion decay 
constant can be substantially lowered from a generic value of a scalar singlet VEV. The 
Yukawa coupling structure in the recent Z12-7 model is studied completely, including the 
needed nonrenormalizable terms toward realistic quark and lepton masses. In this model 
we find an approximate global symmetry and vacuum so that a QCD axion results but its 
decay constant is at the GUT scale. The axion-photon-photon coupling is calculated for 
a realistic vacuum satisfying the quark and lepton mass matrix conditions. It is the first 
time calculation of c a77 in realistic string compactifications: c a77 = I — 1.93 ~ —0.26 
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1. Introduction and Summary on Superstring Axions 

The strong CP problem is, "Why is the QCD vacuum angle #qcd so small at |#qcd| < 
10~ 9 ?" Otherwise, strong QCD interactions violate the CP symmetry and then neutron 
will develop an electric dipole moment of order 10~ 3 |#qcd I x (charge radius of neutron), 
and the present upper bound on the neutron electric dipole moment d n < 0.63 x 10~ 25 e cm 
Q restricts |# QC d| < 10~ 9 . There are a few solutions of the strong CP problem : (i) Set 
Oqcd = at tree level and guarantee that loop effects are sufficiently small, (ii) m u = 
method, and (hi) the Peccei-Quinn (PQ) mechanism. Axion solutions which we discuss 
in this paper belong to Case (iii). The PQ mechanism || introduces an anomalous (in 
the QCD gluon fields) global U(1)pq symmetry. This must be an axial symmetry. Since 
quarks are massive, the global U(1)pq symmetry must be spontaneously broken, generating 



- I - 



a Goldstone boson called axion j| . Currently, the phenomenologically allowable QCD axion 
is a very light axion ||]. A probable initial condition of the axion decay constant allows 
the window, 10 10 GeV < F a < 10 12 GeV. But, with the anthropic principle applied, the 
upper bound can be further open || . 

Axion models in field theory are artificial in the sense that the U(1)pq symmetry is 
given for the sake of the PQ mechanism only. It is desirable if a consistent theory with 
an ultraviolet completion gives a natural candidate for axion. In this regards, we consider 
string models. If a string theory predicts a phenomenologically allowable axion, this may 
be a key prediction of string theory. This is welcome in view of the scarcity of direct 
verifiable methods of 'string' nature. In fact, one attractive feature of string theory is 
the natural appearance of axions from the antisymmetric tensor field Bmn- These super- 
string axions are split into two categories: one is the compactification scheme independent 
one called model-independent axion (MI axion) B^ u |?]] and the other is the compactifi- 
cation dependent one called model-dependent axion (MD axion) Bij ||. However, these 
superstring axions are known to have some problems. Derived from string theory, decay 
constants of these axions are expected to be near the string scale ~ 10 16 GeV ||, outside 
the aforementioned axion window. For MD axions, it is further known that the shift sym- 
metries of MD axions are broken at high energy scales ]!(]] so that they cannot be used 
as axions for rolling the vacuum angles to zero @. However, one should not forget that 
such superpotential of MD axions is a model-dependent statement [11]. Recently, there 
has been attempts to lower the decay constants in string models or in higher dimensional 
models ||, ||, ||, jDj, 0- 0ne approach is a large volume compactification to lower the 



fundamental scale [|T^] . Another approach is using the warped geometry to lower the scale. 
In heterotic flux compactification, MD axions can be localized at a vanishing cycle which 
is warped due to the flux so that we may have a small axion decay constant |fj4[| . Similar 
setup has been discussed in a higher-dimensional model Jlq]. In other contexts, an axion in 
the Kachru-Kallosh-Linde-Trivedi setup has also been discussed in [16]. 

In this paper, we restrict the discussion on superstring axions to heterotic string only, 
but the generic problem of axion mixing is present in any superstring axion models. For 
the MI axion, the decay constant is near the scale ~ 10 16 GeV ||. Some string compactifi- 
cations such as the simple compactifications of Refs. [|l7], [l8|] do not lead to an anomalous 
U(l), in which case the MI axion is harmful ||. Later, it was found that some string com- 
pactifications lead to an anomalous U(l) gauge symmetry |l{J. This anomalous symmetry 
can be gauged owing to the Green-Schwarz mechanism by which the antisymmetric tensor 
field B^ u transforms nonlinearly under the U(l) symmetry so as to cancel the anomaly 
[IS]. This anomalous U(l) gauge symmetry is a subgroup of SO(32) or of EgxEg. The 
Green-Schwarz term |2l|] makes it possible for this anomalous U(l) gauge boson to absorb 
the MI axion and become massive. Below this gauge boson mass scale ~ 10 16 GeV, there 
results a global symmetry U(l) an . This is a kind of 't Hooft mechanism [22|. Under this 
circumstance, if some scalar field carrying U(l) an charges develop VEVs around ~ 10 11 
GeV, then we obtain a harmless very light axion by the PQ mechanism Q. However, 
for this U(l) an , the story is not that simple. Most light fields carry nonvanishing U(l) an 
charges, and if one is forced to give GUT scale VEVs to those singlets for successful Yukawa 
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coupling textures, then the U(l) an breaking scale must be the GUT scale and the resulting 
axion is again harmful. The QCD axion with this U(l) an has been extensively studied 
without any hidden sector confining force in |^] , where the Yukawa coupling textures were 
not used. Moreover, this model is phenomenologically unsatisfactory since sin 2 9y/ <C |. 

Most string compactifications need another confining force in the hidden sector for 
the purpose of introducing supersymmetry breaking. Then, we need an additional global 
symmetry to settle both QCD and hidden sector 9s. As mentioned above, only the U(l) an , 
related to the MI axion, is a good one to consider below the string scale. Except this U(l) an , 
there is no global continuous symmetry resulting from string theory. For an additional 
global symmetry, we are only at a disposal of approximate global symmetries from string 
compactification. It is better for this approximate global symmetry to be broken by a 
sufficiently high dimensional operators in the superpotential so that the symmetry breaking 
superpotential is negligible compared to the axion potential derived from anomalies. This 
idea was examined in a SUGRA field theory model with a discrete Z3 x Z3 symmetry [^4[ . 
But it has not been studied in string compactifications. A discrete symmetry is a good 
guideline to make approximate symmetry violating terms appear at higher orders. If PQ 
symmetry breaking scales are around the intermediate scale, it has been known that the 
PQ symmetry breaking superpotential must be forbidden up to D = 9 terms [25]. But this 
statement is an oversimplified one because higher order terms can involve some scalar fields 
developing small VEVs or even not developing any VEV. So, the PQ symmetry violating 
terms in the superpotential must be checked in mo del- by-mo del bases. Note that for the 
terms breaking the global symmetry to appear at a sufficiently higher order, we need a 
large N, presumably ./V = 12, in the Z^r orbifold compactifications. 

By the way, the MSSMs from superstring need Yukawa couplings beyond cubic terms 
1 26, p?!] . So far, there has not appeared any model where only cubic terms are sufficient to 
give all the needed Yukawa couplings. So it is not unreasonable to require that for realistic 
quark and lepton masses superstring models need nonrenormalizable terms beyond cubic 
terms. 

Thus, to realize a QCD axion in string compactifications, we must satisfy the following 
conditions: 



• One must work in a phenomenologically successful string derived model. 

• One needs an additional confining group beyond QCD, and an approximate global 
symmetry must be introduced. The MD axions cannot be used since the world sheet 
instanton effects violate the shift symmetries of the MD axions. 

• The Yukawa coupling structure must be studied carefully to derive the approximate 
global symmetry. 

• With two axions, the axion mixing effect must be clarified. 

So far, there has not appeared any literature satisfying all the above conditions. In this 
paper, we try to explore the possibility of satisfying all these conditions in the recently 
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proposed string MSSM [p6|1 . We find a vacuum satisfying all of these conditions, but the 
QCD axion decay constant falls in the GUT scale region. 

We emphasize the importance of the last condition which has been overlooked in many 
superstring axion models. It is studied in Refs. 28, 14|. There is the cross theorem on 



axion potential heights and decay constants. It is for the case of two 9s with a complete 
mixing of two axions by the higher potential that the smaller decay constant corresponds 
to the higher height of the axion potentials and the larger decay constant corresponds to 
the lower height of the axion potentials. It is shown for the case of two axions with the 
MI axion and the MD axion corresponding to the breathing mode moduli [29] where the 
couplings are ~ a\{FF + F'F') + a2(FF — F'F'). In fact, the axion mixing occurs when 
both axions couple to the anomaly which give rise to the higher potential from instanton. 
In the example of [29], two anomalies couple to both axions and hence the condition for 
the theorem is satisfied. In such cases, if the hidden sector axion potential is higher than 
the QCD axion potential (as one might guess), then the decay constant of the QCD axion 
is the GUT scale. If we realize that the decay constant of the QCD axion falls around 
~ 10 11 GeV, then it will be observable by axion detection experiments like the CAST of 

cern g. 

After solving the strong CP problem, we can consider a related, the so-called fi problem 
1 31], derivable from the global symmetries of the MSSM. The common origin of a very 
light axion scale and supergravity scale was pointed out early [32]. Phenomenologically, 
the strength of the \x parameter in —^H u H,i is required to be at the electroweak scale. 
The first [i problem is why it is forbidden at the Planck scale. The second \x problem is 
why it is of order the electroweak scale. There have been suggestions that if it is forbidden 
at the GUT scale, it is expected to be generated at the electroweak scale in supergravity 
models [33] and in string models 34]. However, these solutions need some symmetry 
anyway from the Yukawa coupling structure to forbid it at scales below the GUT scale 
1 35]. If many singlet VEVs are required at the GUT scale, then forbidding just H u Hd in 
the renormalizable superpotential as done in |34j is not enough to exclude the /i term at 
the GUT scale. 

In this paper, in addition we try to calculate the axion-photon-photon coupling c a ^^ 
in a realistic superstring model. For the MI axion, a mechanism to find out the global 
U(l) an was explicitly given before [23] but a calculation of c ary in that model has not been 
meaningful because the model is not realistic due to sin 2 6\y <C | and there is no additional 
confining force for supersymmetry breaking. Now we obtained a realistic superstring stan- 
dard model in a Z12-/ construction [26] and here we can calculate c aTy if there exists a 
QCD axion. 

We may introduce two scales (at Mgut and 10 11 GeV) for breaking U(l) an and an 
approximate global symmetry U(l) g iA- The resulting approximate global current must be 
anomalous in the gluon fields and/or hidden confining gauge fields. The axion mixing is 
studied with the hidden sector axion potential. 

In this paper, we pick up an approximate global symmetry U(l) g iA in addition to 
U(l) an . We use a computer program which is based on the Gauss elimination algorithm 
to study all possible U(l) symmetries from any order of Yukawa couplings. A complete 
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analysis has been done for all Yukawa couplings up to D = 9 superpotential terms. Unfor- 
tunately, we could not realize the QCD axion with F a ~ 10 11 GeV but at the GUT scale. 
The axion energy crisis may be resolved by the anthropic principle ||. We calculate the 
phenomenologically important axion-photon-photon coupling constant c a77 , i.e. c a77 — ^ 
or c a77 ~ —0.26. This is the first calculation of c a77 from a realistic string compactification. 

This paper is organized as follows. In Sec. [2], we summarize briefly axion physics 
related to the calculation of the axion decay constant. This involves a discussion on the 
domain wall number of discrete vacua. In Sec. ||, we succinctly present the recent Z12-/ 
orbifold model [26] toward a computer input for the matter fields. In Sec. ||, we summarize 
the Yukawa coupling structure and find the anomalous U(l) an symmetry from EgxEg and 
an approximate anomalous global symmetry U(l) g iA so that two 9s can be settled to zero 
via the PQ mechanism. In Sec. ||, we calculate the axion-photon-photon coupling by 
calculating anomalies. It is compared to the recent CAST experiment bound fl30[ . Sec. || 
is a conclusion. 



2. Axion, Domain Walls and Axion Decay Constant 

The QCD axion a is defined to be a pseudoscalar particle coupling to the gluon anomaly, 

" -G^GT = 4r{GG} (2.1) 



32vr 2 F a ^ ~ F a 

where G^ u is the dual of G^ u , the gluon kinetic term is {\ j \g^)G ^G^ and F a is the axion 
decay constant. It is assumed that there exists the canonical kinetic energy term of a, i.e. 
T^d^ad^a, and there is no potential for the axion except that derivable from Eq. ([2,1|), 
Then, the minimum of the axion potential is at (a) = |3], |3(| which solves the strong CP 
problem cosmologically in the present universe [37]. We do not repeat the explanation of 



the solution of strong CP problem using the axion in this section. For a complete review, 
refer to ||. Here, we discuss some subtlety of determining the axion decay constant due 
to the remaining discrete group after U(1)pq breakdown. 

Because of the quantization of the integral of GG, the axion potential is periodic with 



the periodicity 2isF a . In this form Eq. (^1), the fundamental region of the axion vacua 



is [0,271-^], because 8 = [0,2-7r]. Namely, starting from the vacuum (a) = 0, the next 
vacuum occurs by shifting (a) — > (a) + 2irF a . But the vacuum at (a) = 2irF a may not be 
the same vacuum as the (a) = vacuum, but returns to the (a) = vacuum only after 
the shift a — > a + 2TTNf)wF a . Then there are degenerate vacua whose number is called the 



domain wall number Njjw [38]. The axion embedded in some fields may not return to its 
original value when one shifts a — > a + 2irF a , which is the reason for the appearance of 
degenerate vacua. There is another degeneracy from the scalar field space. Suppose, the 
axion a embedded in the phase of a scalar field as 4> = [( v + p)/V2]e ta l v where <p carries N 
units of PQ charge. The field <p returns to its original value by shifting a — > a + 2irv. But 
the phase factor becomes identity for N distinct points of < a < 2irv. This is related to 
the domain wall number calculation. If we define the PQ charge Q such that <j> carries one 
unit of Q and the anomaly calculation leads to d^Jpq = n{GG}, the phase a couples to 
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I 




Figure 1: Illustration of degenerate vacua for N\ = 4, N2 = 2. Bullets and circles correspond to 
minima of the axion potential, (a) Now is two, the vacua connected by solid and dashed lines, (b) 
The fundamental length (the arrow) in the fundamental region corresponds to the decay constant 
which is F a = V/N DW T. 



the anomaly as ^n{GG}. Since 4> is defined such that it comes to the original value by the 
PQ transformation or the shift of a by a — > a + 2ttv, viz. Q = 1, this interaction defines 
F a = v/n. Here, the domain wall number is Now = n - 

On the other hand, if we are required (by the PQ charges of the other fields) to define 
the PQ charge of 4> as Y which is a positive integer greater than 1, then Jpq = Tn{GG}, 
and we expect the axion decay constant F a = v/nT. But the domain wall number defined 
from vacuum structure is a topological one and still we should obtain Nj^y/ = n. This 
situation is illustrated for N\ = nT = 4 and N2 = T = 2, i.e. with n = 2 in Fig. [I]. 
Certainly, the minima of the axion potential has two distinct sets in Fig. |l|(a), the set of 
bullets and the set of circles, represented by the connections with solid and dashed lines, 
respectively. These two sets agree with our original domain wall number, Now = n = 2. 
This domain wall number is obtained from N\ (the coefficient of anomaly) by dividing it 
with the greatest common divisor of iVi and N%. In Fig. |l](b), we show how our F a is 
related to the original VEV. 

Having specified the domain wall number, it is important to pick up the properly 
normalized axion field toward calculating the axion decay constant. Consider two fields 4>\ 
and 4>2- In the unitary gauge, we can express the Goldstone bosons in the phases of <j)\ and 
4>2 as 

fa = yi + p± ( Mi/vi > ^ 2 = Yi±Pi e *A 2 /v 2 M 

V2 v 2 

whose PQ charges are T\ and T2, respectively. Now we can restrict the phase fields A\ = 
[0, y~Vi] and A2 = [0, j^^]- If we are sitting in one domain, say on the solid line (or on the 
dashed line) of Fig. [l] for reinterpreting the figure with A\ and A2 directions with T\ = 4 
and T2 = 2, A\ and A2 completely define the allowed field space in that domain. Thus, we 
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can represent the fields in one domain as 



<Ai = ^^e^/ 4 ), Al = [0^] (2.3) 
^ = ^fe^ v ^, A 2 = [0,|]. (2.4) 



The relevant VEV for axion is + 2 2 VJ. The axion direction in terms of 'vl's is 

9 A = , 4 ^+ 2A 2 . Suppose tl 

the axion coupling is given by 



^=====. Suppose that the gluon anomaly were just one unit of {GG}. Then, 



4j4l + ^ 2 - r {GG} (2.5) 



V4 2 ^ + 2 2 F 2 2 



With the shifts of A\ — > A\ + 27rVi and A 2 — > A 2 + 271-V2, viz. Fig. |], the field space ends 
at the other domain, i.e. Ai shifted such that 8 changed by 4-7T. So to maintain the fields 
in the same domain with 9 = [0, 2n), we use Ai = 2a{ such that ai — ► ai + 2irFi. Thus, we 
have Fi = Vi/2. Thus, the coupling (|2.5|) can be rewritten as 



1 {GG} = 2(dimensionless axion component in the phase){G(5} (2.6) 



where the overall factor 2 appears as the greatest common divisor (GCD) of T\ and r 2 . 
For one field VEV dominating the other, e.g. the DFSZ model with a large singlet VEV 
(s ) and a small doublet VEV, the dimensionless axion component is in fact ai/(s ) and 
the axion decay constant is F a = (s°)/2. 

If we consider the domain wall number from the gluon anomaly, that should be con- 
sidered also and an appropriate factor must be multiplied to Eq. fl2.6|) . If the anomaly 
coefficient iV(anom.) from the gluon anomaly and iV(phase) from a multitude of phases 
are relatively prime, then the domain wall number is simply iV(anom.). But if they are 
not relatively prime, we must divide it by the greatest common divisor of 7V(anom.) and 
iV(phase). Therefore, a general expression for the axion coupling can be written as 

AW- Dgcd (phase) ^ ^ 

^max 

where -Dgcd is the greatest common divisor of PQ charges of the relevant VEV-carrying 
fields, and the topological number N^w 1S calculated as in Fig. [l|(a) from the coefficient of 
the anomaly and the phase degeneracy, i.e. Now = A r (anom.)/I?GCD(-V(anom.), iV(phase)). 
The normalized axion field is defined as 

° ~ (2.8) 



where Tj is the U(1)pq charge of S{, Vi is the VEV of scalar field (sj), and V max is a typical 
scale of VEVs breaking the global symmetry. In many cases it is the maximum value among 
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(sj). For example, for the case V ma _ x ^>Vi{i^ max) and also for V max = V^(for any i). Then 
the axion decay constant is 

p a = Vmax . (2.9) 

Note that if Nqw ■ Dqqd were large, then F a can be lowered significantly compared to the 
PQ symmetry breaking scale. 

The discussion of this section will be used in Sec. |5] in estimating the decay constant 
of the QCD axion from superstring. 



3. Z 12 _ 7 Orbifold Model 

In this section, we briefly summarize the key points of a realistic string-derived minimal su- 
persymmetric standard model (MSSM) from a Z12-/ orbifold compactification ^6|. Here, 
the low energy effective theory contains the MSSM. We discuss vacuum configuration of the 
MSSM singlets for realistic low energy phenomenology, especially toward the absence of ex- 
otic matter spectra and general formulae for realistic Yukawa couplings. These formulae are 
used in the program to generate all Yukawa couplings. We pick up the anomalous gauged 
U(l), U(l) an , and approximate global U(l)s which appear in general in lower dimensional 
superpotential terms. For gauged U(l)s, the gauge U(l) directions can be assigned in the 
EgxEg group space but for the approximate global U(l)s it cannot be represented in that 
way. The only method for approximate global U(l)s is to list the U(l) charges of the 
matter fields. 



3.1 Flipped SU(5) matter spectrum and Yukawa couplings 

In the orbifold compactification in heterotic string theories, the action of twisting is rep- 
resented as a shift vecdtor V^I = 1,2, ... , 16) and Wilson Lines a\ (i = 1,2, ... ,6; I = 
1, 2, . . . , 16). In the Z12-/ orbifold model discussed in 26], we choose SU(3)xSO(9) lattice 
for toroidal compactification and 

^=(Hiii^n°)(^n ) (3.1) 
a 3 = a 4 = (0 5 § ) (0 § 5 ) (3.2) 

0>l = 02 = 05 = &6 = 0. 

Then, below the compactification scale one obtains the following unbroken gauge group 

[SU(5) x U(l) x x U(lf] x [£17(2) x 50(10) x U(l) 2 }'. (3.3) 

The spectra for the unbroken subgroup SU(5)xU(l)x is in fact the flipped SU(5) model 
]26[| . The matter representations under the nonabelian gauge groups SU(5), SU(2)' and 
SO(10)' are simply determined by embedding the momenta P of EsxEg group space into 
the weight space of the corresponding subgroup. The U(l) charges for matter fields are 
determined by 

q i = Z l -(P + k(V + m f a 3 )), (3.4) 
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where i = X, 1, 2, 3, 4, 5 denotes each U(l) gauge group and k corresponds to the orbifold 
twist number of k-th twisted sector and m/ is the Wilson line twist number applicable only 
to the second two-torus, viz. ( |3.2|) . Zi are the following EgxEg weights 



Zx = (2,2,2,2,2;0 3 )(0 8 )' (3.5) 

Z 2 = (0 5 ;1,0,0)(0 8 )' (3.6) 

Z 3 = (0 5 ;0,1,0)(0 8 )' (3.7) 

Z 4 = (0 5 ;0,0,1)(0 8 )' (3.8) 

Z 5 = (0 8 )(1,1;0;0 5 )' (3.9) 

Z 6 = (0 8 )(0,0;1;0 5 )'. (3.10) 



We will rearrange these U(l)s when we consider the gauge anomalies of the model com- 
pletely. 

The 4D Lorentz symmetry representation is solely determined by the right-mover 
oscillators. Since we have supersymmetry in the low energy limit, we consider either 
the fermion spectrum or the boson spectrum. Consider the fermion spectrum. It is 
determined by the Ramond sector vacuum of the right movers which is SO (8) spinor 
s = (sq,s) = {±^,±^,±^,±^} with an even number of minus signs. The 4D chirality 
X is the first component of s, and we call X = \ the- right-handed field and X = ~\ the 
left-handed field. 

Matter spectrum must satisfy the mass-shell conditions : 

j 

3 

where j runs over {1,2,3,1,2,3} and 4>i — k<f>i and (fc = —k<f>i where = denotes one plus 
the maximum integer smaller than the original real number. Here and are the 
oscillator numbers, and c and c are the zero point energy of the right mover and left mover, 
respectively. 

By the modular invariance, the matter spectrum satisfying the above mass shell con- 
ditions can be further projected out and have multiple copies of themselves. This can be 
summarized by one-loop partition function. The correct formulae are reviewed in | |40| ], 




(3.13) 



where N is the order N in the Z^r orbifold, N\y is the order of the Wilson line, 3 in our 
case, and 

% = " N J R ^i ~ ^ V f ~ ^ + (P + kVf) -V-(s + k4> s ) ■ cf> s , (3.14) 

j 
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(3.11) 
(3.12) 



Name 


SU(5) u{1)x 
xSU{2)' 


U(lf 


P 


Name 


SU(5) u(1)x 
xSU(2)' 


[/(!)* 


P 


Untwisted 


Twisted 6 


U3 5 


(5 3 ,1) 


(6,6,6,0,0) 


1 


T6 5 
T6 5 


(5 3 ,1) 

(5-3,1) 


(6,0,0,0,0) 
(-6, 0, 0, 0, 0) 


2 
2 




(10-1,1) 


(6,-6,-6,0,0) 


1 


^6to 
T610 


(10-i,l) 
(10i,l) 


(6,0,0,0,0) 
(-6,0,0,0,0) 


4 
3 


U3i 


(1-5,1) 


(6,6,6,0,0) 


1 


T6i/,T6u 




(±6,0,0,0,0) 


2,2 


U2- 5 


(52,1) 


(-12, 0, 0, 0, ) 


1 


hi, hi 


(lo,l) 


(0, ±6, ±6, 0,0) 


4,2 




(lo,l) 


(0,0,0,24,0) 


1 


h 2 ,h 2 


(lo,l) 


(0, ±6, ±6, 0,0) 


2,3 


Ul ro 


(10-1,1) 


( 6, 6, 6, , ) 


1 


h3,h~3 


(lo,l) 


(0, ±6, ±6, 0,0) 


2,4 


Ul 5 


(5 3 ,1) 


(6,-6,-6,0,0 ) 


1 


/14 , /14 


(lo,l) 


(0, ±6, ±6, 0,0) 


3,2 


Uh 


(1-5,1) 


(6,-6, -6,0,0) 


1 










Twisted 1 


Twisted 7 




(5_i,l) 


(-7,6,0,4,0) 


1 


T7 M 


(5_i,l) 


(-7,0,-6,4,0) 


1 


Tl- 5B 


(5_i,l) 


(5,-6,0,4,0) 


1 


T7~ 5B 


(5_i,l) 


(5,0,6,4,0) 


1 


Tl 5 


(5i , 1) 


(-1,0,6,4,0) 


1 


T7 5 


(5x,l) 


(-1,-6,0,4,0) 


1 


T1 1A 


(1_|,1) 


(-7,-6,0,4,0) 


1 


T7 1A 


(1-|,1) 


(-7, 0, 6, 4, ) 


1 


T1\B,C 


(1-|,1) 


(5,6,0,4,0) 


1 


T7\b,c 


(1_|,1) 


(5,0,-6,4,0) 


1 


T\ XD 


(1|,1) 


(11,0,6,4,0) 


1 


T7 1D 


(1|,1) 


(11,-6,0,4,0) 


1 


TliE,F,G 


(1|,1) 


(-1,0,-6,4,0) 


1 


T7\e,F,G 


(1|,1) 


(-1,6,0,4,0) 


1 


Si 


(1_|,2) 


(5,2,4,-8,-4) 


1 


h 


(1_|,2) 


(5,-4,-2,-8,-4) 


1 




(1_|,1) 


(5,2,4,4,8) 


1 


T7t B 


(1_|,1) 


(5,-4,-2,4,8) 


1 




(1| 2 ,2) 


(-1,-4,-2,-8,-4) 


1 


$4 


(1|,2) 


(-1,2,4,-8,-4) 


1 


Tito 


(1|,1) 


(-1,-4,-2,4,8) 


1 


T7t D 


(1|,1) 


(-1,2,4,4,8) 


1 


Tl 5 " 


(5|,1) 


(-1,4,2,4,-8) 


1 


T7 5 - 


(5,,1) 


(-1,-2,-4,4,-8 ) 


1 


T^a 


(1_|,1) 


(5,-2,8,4,-8) 


1 


T7- lA 


(1_|,1) 


(5,-8,2,4,-8) 


1 


TKb 


(1_|,1) 


(-7,-2,-4,4,-8) 


1 


T7- lB 


(1-|,1) 


(-7,4,2,4,-8) 


1 


Tl w 




(-1,-8,2,4,-8) 


1 


T7 W 




(-1,-2,8,4,-8) 


1 



Table 1: Spectrum of the Kim-Kyae Zi 2 _/ model: 1. U sector and T6, XI, T7 sectors. 



where (f>i = 4> s iSgn((f>i) and Vf = (V + m/a). Here, x{9 k ,6 l ) is the degeneracy factor 
summarized in p6]| . Vk is interpreted as the multiplicity of the spectrum. Thus, if Vk = 0, 
then the field is projected out by the modular invariance. Multiple Vk implies that the 
same matter spectra occur at different orbifold fixed points. We summarize the resultant 
matter spectrum of this model in Tabel [l] and 0. 

We summarize the left-handed matter field spectra, 1 by classifying them under the 
1 We show only the left-handed spectra. The right-handed spectra are just their PCT conjugates. 
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Name 


^(5)c/(i )x 
x S*7(2)' 




P 


Name 


SU(5) V ( 1)x 
xSU{2)' 


[/(I) 5 


P 


Twisted 2 


Twisted 4 




(5s, 1) 


(-2,0,0,8,0) 


1 


T4 5 


(5-2,1) 


(-4,0,0,-8,0) 


3 


T2i 


(1-5,1) 


(-2,0,0,8,0) 


1 


T4g 


(52,1) 


(-4,0,0,-8,0) 


2 


^1,2 


(lo,l) 


(4,6,-6,8,0) 


1 


„o 

1,2,3 


(lo,l) 


(8,0,0,-8,0) 


2 


u 3,4 


(lo,l) 


(4,-6,6,8,0) 


1 


„0 
s 4 


(lo,l) 


(8,0,0,16,0) 


2 


r o 


(lo,l) 


(-8,6,6,8,0) 


1 


„0 
s 5 


(lo,l) 


(-4,12,0,-8,0) 


2 


c 6 ° 


(lo,l) 


(-8,-6,-6,8,0) 


1 


r\ 

„U 

S 6 


(lo,l) 


(-4,-12,0,-8,0) 


2 




(lo,l) 


(4,-6,6,-16,0) 


1 


n 

s 7 

n 

s 8 


(lo,l) 
(lo,l) 


(-4,0,12,-8,0 ) 
(-4,0,-12,-8,0 ) 


2 
2 


T2 OT6 


(lo,l) 


(4,-2,2,-4,-2) 


1 


df 


(lo,2) 


(8,-4,4,4,-4) 


2 




(lo,2) 


(4,2,-2,-4,-4) 


1 


4 


(lo,2) 


(-4,8,4,4,-4) 


2 


^1,2 


(lo,l) 


(4,2,-2,8,8) 


1 


4 

4 
4 


(lo,2) 
(lo,l) 
(lo,l) 
(lo,l) 


(-4,-4,-8,4,-4) 
(8,-4,4,-8,8) 
(-4,8,4,-8,8) 

(-4,-4,-8,-8,8) 


2 
3 
2 
2 




(lo,l) 


(4,2,-2,8,-4) 


1 


di 


(lo,2) 


(8,4,-4,4,4) 


2 




(lo,2) 


(4,-2,2,-4,4) 


1 


d 2 


(lo,2) 


(-4,-8,-4,4,4) 


2 


^1,2 


(lo,l) 


(4,-2,2,8,-8) 


1 


d 3 


(lo,2) 


(-4,4,8,4,4) 


2 




(lo,l) 


(4,2,-2,-16,8) 


1 


s i 

s 2 
s 3 


(Id) 
(Id) 
(lo,l) 


(8,4,-4,-8,-8) 
(.4,-8,-4,-8,-8) 
(-4,4,8,-8,-8) 


3 
2 
2 



Table 2: Spectrum of the Kim-Kyae Z\2 model, continued: 2. T2 and T4 sectors. Here, 016 and 
O10 mean 16 representation and 10 representation of SO(10)', respectively. 



flipped-SU(5) and SO(10)', 

U : (1-5 + 53+10-1)%, (5 2 )i/ 2 , (I-5 + 53+10-1)%, (1 )%, (3.15) 
from the untwisted sector, and 



n 


T0_i + {2(1_5 + 1 5 + 5_3 + 5 3 ) + 3(T0_i + 10i)} + 22{1 }, 


(3.16) 


T 2 


1 5 + 5 3 + 11{1 } + 4D + O10 + OT6, 


(3.17) 


T 4 


5-2 + 2(5_ 2 + 5 2 ) + 30{l o } + 12D, 


(3.18) 


Ti 


2(5 i) + 2(5, i) +6(1 5) + 6(l + 5) + (Dl +5 ) + (Dl s), 

2 '2 2 '2 '2 2 


(3.19) 


T 7 


2(5 + i ) + 2(5_i) + 6(l + s) + 6(l_s) + (Dl +5 ) + (Dl_t), 

2 2 2 2 2 2 


(3.20) 



where D denotes a doublet of SU(2)', and O10 and 016 mean 10 and 16 of SO(10)', 
respectively. From the hidden sector, there are twenty SU(2)' doublets and one 16 and 
one 10' of SO(10)'. Note that matter fields from Tl sector and T7 sector have exotic 
representations which are not present in the MSSM. We name them by G-exotics for 
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(anti)-fundamental representations of SU(5) and by E-exotics for singlets of SU(5) but 
with nonvanishing U(l)x so that they have exotic U(l) em charge. 

The superpotential terms are obtained by examining vertex operators satisfying the 
Z12— i orbifold conditions [41, |40| ]. It can be summarized by the following selection rules: 

• H-momentum conservation with (f> s = 

^fii(z) = -lmod 12, ^i? 2 (2) = lmod 3, ^ R 3 (z) = 1 mod 12, (3.21) 

Z z z 

where z{= A, B,C, . . .) denotes the index of states participating in a vertex opeartor. 

• Space group selection rules: 

^fc(z)=0mod 12, (3.22) 

z 

[kmf] (z) = mod 3. (3.23) 



Neglecting oscillator numbers, the i/-momenta for Zi2_/ twist are 
U x : (-1 0), U 2 : (0 1 0), U 3 : (0 1), 

r i: (^&&), T 2 :(^|i), T 3 :(i0f) (3.24) 

T 4 : (if I l)> T 5 : (n TJ if )> T 6 : (tt l) 

Now we present superpotential terms relevant for the low energy phenomenology. Such 
terms include the mass terms for the exotics, the \x terms and the Yukawa couplings. We 
summarize them in Tables ^ and 

3.2 Hidden sector SO (10)' and SU(2)' 

There are two hidden sector nonabelian groups, SO(10)' and SU(2)'. Inspecting SU(2)' 
matter spectrum, there are twenty SU(2)' doublets: Si, 82, £3, £4, Dj , Z?{" , D^, , d\, 
rff, d^ , d^ , d% , and d%. Inspecting the SO(10)' matter representation, we find only a single 
SO(10)' spinor 16' and a single SO(10)' vector 10'. Since there are four nonabelian gauge 
groups in total after the flipped SU(5) breaking, it seems that we need to consider four 6 
parameters, those of SU(3) C , SU(2)vy, SO(10)' and SU(2)'. In fact, among those, we need 
to consider only SU(3) C and SO(10)' because S\J(2)w is broken at the electroweak scale 
and we assume that SU(2)' is broken. For SO(10)', here we just assume that a subgroup 
of SO(10)' confines at an intermediate scale since the hidden sector dynamics is not well 
understood yet at present. 



4. Approximate U(1)pq Symmetry of Z 12 / Model 

In this section, we find U(1)pq symmetries allowed from the superpotential terms in the 
Z12-1 model. To all orders, only gauge and discrete symmetries can remain valid since 
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Mass terms for vector-like representations 



T6i T6 f0 1 
T65T65 > 
{ T6i/T6u J 



[ r i ^r , ^n~ to+ TO+ to - 
1 o 5 o 6 o 3 j i QT ^i z oJS i z ow , 



'4i 



a u „0 
6 6 2 , 



4( h 3 h i + s lhih 3 ), 
s u s° 2 (h 3 h 1 + ^3), 
[C 6 °^°(C 6 % + C 5 %) 



;U terms 



T4 5 C/2g 



x < 



„0„0 



„tt „0 „0 

6 6 2 6 2 , 



52^4(^-3^-1 + h\h 3 



s u s° 2 s^(h 3 h 1 + ^3), + hjii), D+D 2 s 4 (h 2 h 3 + h 3 /t 2 ), 

c s4(44 h i + + + dtd^Jn), 

C°s 2 s° 2 slh 3 , Cls^slslhs, Cls14dih 3 , 
Cls^sl^Cls^sfhs, Clsld+d^h 



Mass terms for G-exotics (Ti5,Tl 5 , T7s,T7 5 , TTg^, T1§ B , T7$ A , T7§ b ) 



r7 5 -ri 5A x 

T7^T7 M x 
T7^Tl- 5B x 



sjh 3 

C 3 C<5 (44_ s l + ^^s" 5 ! ) 
s+(l + fr 3 /ii+/nfr 3 ) 



^3^5 ( S 5 S 8 S 1~ + S t S 3 S l ) 



n 5 -T75 A x 

T7 5 Tl- 5A x 
T7 5 T7 M x 
T7 5 Tl- 5B x 

T7 5 T7 EB x 



sJil + h^s + hsh!) 
sl{l + h 3 hi + hji 3 ) 

+ C§s?sg + D+D 2 h 2 ) 



Cis° 5 (s 8 h 3 + s° 5 h 3 



( s^s^/ii/ii + s^hihi 
3 I +44hihi 



Mass terms for E-exotics 



Tl+ c T7 1A x 
Tl 1E T7 1B x 
Tl 1B Tl 1D x 
T7 1B T7 1D x 
T7 1D T7+ B x 



d 1 (1 + h 3 hi + hih 3 ) 



(CgSgSg + C'gSg" S 2 

(Ce44 + C6 S 2 s 3 



T1 1A T7+ D X S^(h 3 h 2 + h 2 h 3 ) 



T\ 1E T7~ B x 
Tl 1B Tl 1E x 
Tl+ D T7 1A x 
T7 1A T7 1G x 
T1 1A T7+ C x 

T7 lB T7+ D x 



s+(l + h 3 h 1 + h 1 h 3 ) 

(Si 7*3^2 + Sl ^2^3) 

h/t(C 2 Cj + C 2 C 3 ) 
d\[h 2 h 2 



I C l 2' s 7' s 2 + C6 S 1 
[ +C 2 SgS 3 + C 2 



2 S 2 S 3 



Table 3: Superpotential terms relevant for phenomenology: f . Mass terms for vector-like repre- 
sentations (^i-like terms), G-exotic fields and E-exotic fields, and the \i terms for H u and Hd- We 
represent superpotential terms up to dimension 5 only for mass terms of E-exotic fields since there 
are too many terms from dimension 6. 



string theory a priori does not impose any continuous global symmetry. However, at a cer- 
tain order of the effective Lagrangian, the theory can have approximate global symmetries. 
Such approximate symmetries are not given from the first principle, and thus we have to 
enumerate allowed symmetries by reading all superpotential terms. To do this, we have 
made a computer program which enumerate all the allowed potential terms and continuous 
abelian symmetries at a given order. 

If we assign U(l) charge qj for each field <£j, a superpotential term Wi oc f7 ^ ( <I>j ) Ai J has 
U(l) charge ^ ■ KjQj, where A^- represents the number of occurrences of the field $j in the 
i-th superpotential term. Thus, the algorithm for finding out U(l) charges is basically the 
same as finding eigenvectors with zero eigenvalue of the matrix X T X. It is remarkable that 
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Yukawa Couplings 


Generation 


Matter 


i 


17%, UU, UU 




o 

z 


U3 W 


, C/3 5 , ?73i 




Q 

o 






r 2i 




(ij) 


rojiOj-Bn 


(ij) 


10^5^5// 


(3,3) 


T6 To T6 To U2 E x {l,/^,^} 


(3,3) 


T6 


1 - T2 5 r4 5 x {1,/^iAM 


(3,2) 


TGioUS^-sh! 


(3,2) 




T6 f0 ?735r45C 6 


(3,1) 


T6 To Ul ro U2sh 3 


(3,1) 




TGroUhT^Cl 


(2,2) 


U3^u3^u2- b h 1 h 1 


(2,3) 




U3 m T2 5 T4 5 h 1 


(2,1) 


U3 f0 Uli U2- 5 x {l.As/ii./nha} 


(2,2) 




U3 fG U3 b TA b C^h 1 


(1,1) 


C/lfof/lro^g/ta/Js 


(2,1) 




C/3 f0 C/l 5 T45C 5 /ii 






(1,3) 




[/l 1 - T2 5 r45^3 






(1,2) 




Ul f0 U3 5 T4 5 Cfh 3 






(1,1) 




Ul ro UhTA 5 Cfh 3 


(ij) 








(3,3) 


T2 5 T2 l U2- 5 s%s° 2 








(3,2) 


T2 5 U3 1 U2 E s% 








(3,1) 


T2 5 Ul 1 U2- 5 s° 2 h 3 








(2,3) 


U3 b T2 1 U2is° 2 h l 








(2,2) 


C/3 5 C/3i?725/ii/ii 








(2,1) 


E/3 5 E/litf2 g x {1,^1, Ma} 








(1,3) 


[/l 5 T2 1 ;725 S 0/ l 3 








(1,2) 


C/l 5 C/3i?725 x {l,/i 3 ^i,/ii^3} 








(1,1) 


I71 5 ?71i?72s/i3/i3 









Table 4: Super potentential terms relevant to the low energy phenomenology continued: 2. 
Yukawa couplings of MSSM matter fields. 



we can find the eigenvectors exactly (without any numerical error) since the eigenvectors 
can always be represented by integer-valued vector due to the fact that the matrix A T A has 
only integer-valued components and we are only interested in the eigenvectors with zero- 
eigenvalue. Eigenvectors with zero eigenvalue can be obtained by the Gauss elimination 
method: By multiplying 



A 



or 



V 



(4-1) 



■I 



onto A A, the eigenvector solutions are not changed. Using such operations, one can reduce 



-14- 



the matrix up to the following upper triangular form 2 

/ on ai2 ■■■ clin \ 
'•■ ••• : 



(4.2) 



\ ••• a NN / 

where N is the total number of field species, the diagonal element an is either 1 or 0, 



dij = when 



and ciji = when an = 1. Denote nonzero element of ay by b 



r runs over k indices of nonzero diagonal elements and a runs over N 
diagonal elements. Then, the eigenvectors are 



Ur\q a } 



{(-brl) 
K-br2) 



|100 ... 0} 
|010 ... 0} 



ra where 
k indices of zero 



(4.3) 
(4.4) 

(4.5) 
(4.6) 



{Hv(Ar-fc))|o...i}, 

where q r are the eigenvector components of the r-th element with nonzero diagonal element 
for the above matrix, and q a are those of zero diagonal element. 

We have found all the superpotential terms up to dimension 9. From these, we have 
found U(l) symmetries using the algorithm described above. Certainly, six of them are 
U(l) gauge symmetries. Obtaining all the gauge anomalies, we rediagonalize the gauge 
symmetries according to the anomalies: 

' Qx = -Zi 
Qi = Z 2 + 6Z 3 

Q 2 = -Z 2 + 6Z 4 (4.7) 
Q3 = z 5 
Q4 = 2Z 2 + SZg 

Qan = -6Z 2 + Z 3 - Z 4 + 4Z 6 (4.8) 

Five U(l)s of (p~7|), U(l)x,U(l)i, • • • , and U(l)4, do not carry any gauge and grav- 
itational anomalies while the U(l) an of ( [4.8|) is anomalous. Therefore, this model has 
an anomalous gauge U(l) which can be consistent only with Green- Schwarz mechanism 



1 21, 19 1 . Note that the anomaly check is very nontrivial, and hence this gaurantees the 
matter spectrum shown in Table and ^ is correct. 

Up to dimension 7, we find one global U(l) symmetry which does not belong to the 
above gauge symmetries. This accidental global symmetry must be broken at still higher 
order superpotential terms. We discuss it in Subsec. |4^ 

The anomalous gauge U(l) induces the Fayet-Illiopoulos term. D-flat condition for the 
anomalous U(l) an is required to be 



(D 



(an)\ 



2.9 



192vr ; 



: TrX an + ^X 

anVVr V 



I 0(1 



0. 



2 Note that we cannot make the matrix be the identity since the matrix is not invertible. 
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Some singlet fields get VEVs to satisfy the above equation, breaking some U(l) symmetries. 
But we will not focus on the details of this mechanism. Below, we try to find a vacuum 
(or vacua) realizing observable QCD axion. For this objective, certainly we will choose 
some VEVs of singlet scalar fields. In doing so, we will restrict such that all the known low 
energy phenomena such as fermion masses are reproduced. So, we will consider Q em = |, — i 
quarks and charged leptons. Neutrino masses are not considered here since the mass matrix 
is too gigantic because of the appearance of numerous singlets. 



fields 




fields 




fields 




fields 




fields 


Qan 


U3 5 


-3 


1U 


-3 


U3 1 


-3 


U2 E 


6 


U1 T0 

1U 


-3 


Ul 5 


-3 


Ul x 


-3 


TIm 


4 


TI5B 


-3 


Tlu. 


3 


T1 1B 


-2 


Tl lc 


-2 


T1 1D 


-6 


TllE 


1 


T1 1F 


1 


T\ XG 


1 


Si 


-4 


h 


-1 


Tl tn 


3 


T\7 




-2 


T1 1A 


-6 


T1 1B 

Lid 


1 


Tl ic 


-3 




1 


T2i 


1 


T2 OTE 


-3 


1 z O10 


-3 




-1 


cl 


-1 




-3 


c\ 


-3 




4 


r<o 

°6 


4 




-3 


DT 


-1 


D 2 


-1 


cr 


-5 


C 2 " 


-5 


C3 


1 


Dt 


-3 


Dt 


-3 


ct 


1 


ct 


1 


T4 5 (3) 


2 


T4g(2) 


2 


«°(2) 


-4 


4(2) 


-4 


■m 


-4 


4(2) 


-4 


4(2) 


3 


■8(2) 


1 


4(2) 


1 


4(2) 


3 


4(2) 


-6 


4(3) 


-2 


4(2) 


1 


4(2) 


5 


4(2) 


1 


4(2) 


5 


4(2) 


-2 


*r(3) 


-6 


4(2) 


3 


^ (2) 


-1 


4(2) 


3 


*3~(2) 


-1 


T6 5 (2) 


-3 




-3 


T6 U (2) 


-3 


T6g(2) 


3 


T6io(3) 


3 


T6u(2) 


3 


T7~ 5A 


4 




-3 


T7ia 


3 


T7 1B 


-2 


T7 1C 


-2 


T7 1D 


-6 




1 




1 


T7i G 


1 


S3 


-4 




-1 


J 'id 


3 


T7 5 


-2 




-6 




1 


T7 W 


-3 















Table 5: Non-vanishing U(l) an charges. Fields ^{V) denotes the multiplicity V in case V > 1. 



Low energy fields and the anomalous gauge U(l) an charges are listed in Table ||. D, d, 
and 5 fields are SU(2)' doublets. One can easily check that the U(l) an -SU(5)-SU(5), 
U(l) an -SU(2)'-SU(2) / , and U(l) an -SO(10)'-SO(10) / anomalies are universal, which is re- 
quired by Green- Schwarz mechanism. 

We must assign some singlets GUT scale VEVs for successful Yukawa couplings. First, 
the singlet fields combinations which appear in the Yukawa couplings must have GUT scale 
VEVs generically. Since the exotic fields which appear in T\ and T7 sectors must have 
the GUT scale mass not to spoil the gauge coupling unification, the singlet fields in those 
mass term must be of the order the GUT scale also. 

4.1 Exotics 

There are two kinds of exotics, G-exotics and E-exotics. They are named according to 
SU(5): G-exotics are SU(5)x (anti-)quintets and E-exotics are SU(5) singlets (with non- 
vanishing X charges). 
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We will show shortly in Subsec that should be small to explain the electron mass 
successfully. Under this circumstance, we try to remove exotics at the GUT scale. 

4.1.1 G-exotics 

We have 4 pairs of G-exotics from Tl and T7 sectors: 

5 1/2 :T1 5 ,T1+,T75,T7 5 - 
5-1/2 : T1^,T1^ B ,T7^ A ,T7^ B 

In terms of SM gauge group SU(3)c x SU(2)i x U(1)y, they transform as 5i/ 2 ~^ (3, l)i + 
(1, 2)q and 5_ 1 / 2 — ► (3, + (1, 2)o- In Table |3[ we show the mass terms for G-exotics 

6 

up to dimension 7. Some mass terms at dimension 6 and 7 are omitted if there are lower 
dimensional terms of the same type since they are too many. All those mass terms which 
do not appear in the table have terms at dimension 8 except for T1^T1^ B and T1^T7^ B . 

4.1.2 E-exotics 

There are sixteen E-exotics. It is not wieldy to calculate the determinant of the E-exotics 
mass matrix. Therefore, we first find out approximate U(l) directions from other Yukawa 
couplings and then look for the E-exotics mass matrix. 

4.2 Vectorlike pairs 

Another consideration is the mass terms for the following vector-like pairs, 

(T65,r65),(r6 1 o,T6 T o),(r6 1/ ,r6ij), (T4 5 ,T4g). (4.9) 

The vector-like pairs in T6-sector of ( [4.9| ) are exactly vector-like under all gauge symmetries. 
Also, they have exactly the same i?-momenta, which means that all the vector-like pairs 
(T65, T6g), (T610, T6fo), (T6ij, T6u) satisfy the same selection rules for the superpotential. 

To break the flipped SU(5) down to the SM gauge group, vectorlike pairs 10 + 10 of 
T6 must develop VEVs. Here, note that the multiplicities are different: 'P(IO) = 4 and 
"P(IO) = 3. One 10 is not matched and it becomes the third family member called 10< 
of the SM. Thus, the coupling allowed by T6T6 is restricted to three vectorlike pairs of 
10 + 10, and 10^ carries an independent phase. At the lowest order, the resultant mass 
terms for three vectorlike pairs are presented in Table ^. Those two terms involve the 
hidden sector matter fields T2~t— and T2^ in . T2~t— cannot have a VEV at the GUT scale 

do IVJ.U do 

since it leads to the D-term SUSY breaking at that scale. A VEV of T2q 10 has no problem 
with SUSY breaking, but here we do not assume about the symmetry breaking pattern 
in the hidden sector. At dimension 9, which is next to the lowest order, there are many 
terms inducing mass terms for the vectorlike pairs. In this paper, we do not show all the 
dimension 9 terms since there are too many of them. Since we will argue that C5 and CS 



must have a VEV near the GUT scale as shown in Subsec. [4.4| , we show the mass terms 
involving C5 and C§: 

T6i T6fo 1 f + 01, 1 + OS 1 

T61/T61J I 3 3 2 3 > 
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Therefore, if a combination of (C3, s 1 , or sj, or S3) has a O(Mqut) scale VEV, the 
vectorlike pairs can have O(Mgtjt) scale masses. 

4.3 Approximate global symmetry U(l) g iA 

Before proceeding to discuss the global symmetry we have found and used it as a PQ 
symmetry for QCD axion, we first discuss the difficulty of U(1)pq with the intermediate 
PQ symmetry breaking scale around 10 10 ~ 12 GeV in ^12-/ orbifold flipped SU(5) model. 

The following fields are non-singlet under SM gauge group. MSSM matter and Higgs 
fields : 

Ul m f/l 5 , Uh, 

17%, U3 5 , U3 U (4.11) 
T6 io> T2 5, T2 U 

Th,U2- 5 , (4.12) 

Vector-like pairs of SM matter-like representation : 

(T4 5 ,t4 5 ), 

(T6 10 ,T%), 

(T6 5 ,T6g), 

(T6u,T6u), 

and vector-like pairs of G-exotics and E-exotics from the 1st and the 7th twisted sectors. 

One pair of T610 and T6jq must have GUT scale VEV to break the flipped SU(5) 
gauge group to the SM ones and then T6s,T6g and TQuT6\j must also have the mass of 
O(Mqut) not to spoil gauge coupling unification. On the other hand, 74s and T4§ do 
not need to have the GUT scale mass since they form a complete SO(10) multiplet. Note 
that the nonzero PQ charges of T610 an d T6jq do not necessarily mean PQ symmetry 
is broken at the VEV scale of those fields because the Goldstone boson of the broken 
symmetry is eaten as a longitudinal degree of freedom of the massive gauge boson and one 
linear combination of the original global symmetry and one of SU(5) x U(l)x can remain 
as a surviving low energy global symmetry, which we called the 't Hooft mechanism p2| . 
However, to have such linear combination exist, T610 an d T6jq must have a vector-like pair 
of charges under U(1)pq symmetry. Furthermore, the other fields which become massive 
at the GUT scale must have vector-like charges since we want to break U(1)pq at much 
lower scale. 

We may try U(1)pq symmetry to be flavor-dependent. However, then the flavor- 
dependence is severely restricted if we consider the PQ breaking symmetry scale around 
0(1O 12 GeV) since the maximum hierarchy we have in the low energy Yukawa coupling is 
10~ 6 for the Q e m= § quarks. 
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Firstly, let us consider the flavor-independent PQ symmetry case. Denote U(1)pq 
charges of SM non-singlet fields by 

Qio = Q(U1 W ) = Q(U3 W ) = Q(T6 W ), 
Q 5 = Q(Uh) = Q(U3 5 ) = Q(T2 5 ), 

Q 1 = Q(U1 1 )=Q(U3 1 ) = Q(T2 1 ), (4.13) 
Q 5H = Q(T4 5 ), 
Q- 5H = Q(U2- 5 ) 

The cubic Yukawa couplings 10^10^5//, 10j5j5// and 5jlj5# are required to be invariant 
under the symmetry we try to find, and we obtain the following relations, 

QbH = —Q\o — Q5, 

Q 1 = 2Q T5 -Q 5 . (4.14) 

In addition to the SM Yukawa coupling, this model has another cubic coupling involving 
SM non-singlets: 

T65T45T21 + T6 5 ?72gT6i/ + T2 5 T4^T6 U - (4.15) 

From the above couplings, we further get 

Q(T4 5 ) = Q- 5H = -Q w , 
Q(T6 5 ) = -Q(T6g) = Q 5 , 

Q{T6 U ) = -Q(T6u) = Q 1 = 2Q m - Q 5 . (4.16) 

Since this determines all the U(l)pQcharges for SU(5) nonsinglets, we can calculate the 
U(1)pq-SU(5)-SU(5) anomaly in terms of Qjq and Q5. One can easily check that it is just 
zero, and thus we cannot have any approximate flavor-independent PQ symmetry which 
has the U(1) PQ -QCD-QCD anomaly. 

This leads us to the flavor-dependent case. However, the situation does not improve 
very much in this case either. Since the down type Yukawa coupling 10^10^5^ does 
not reveal large hierarchy among the families, we have to assign the flavor-independent 
symmetry for the field which take part in these Yukawa couplings. Thus, it still holds that 

Qio = Q(Ulw) = Q(UZw) = Q(T&ro). (4.17) 

Additionally, trilinear couplings give a strong restriction in choosing a PQ symmetry. If 
a cubic combination of SM nonsinglet field is not invariant and that combination appears 
with other singlet fields in the superpotential, one can assign the compensating U(l) charges 
for the involved singlet fields so that the overall term is still invariant. However, for the 
trilinear couplings of the SM fields which does not involve any singlet fields, this cannot 
be done. Therefore, such terms must be required to be invariant first. In the SM Yukawa 
coupling combination, we have 

T6 I qT6 I ( J [/25 + C/3 T oC/l T oC/25 + T^T2^TA b 

+U3 b Ul 1 U2i + Ul b U3 1 U2i. (4.18) 
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Another consideration we should take care of is the overall mass hierarchy of the SM 
matter fields. From the observed masses of quarks and leptons, we conclude that the 
determinant of up- type quark mass matrix and lepton mass matrix must have at most 
one factor of M ^ ® '- . This requires at least one of 22-component and 33-component of the 
up-type quark Yukawa couplings must have no U(1)pq charge. In fact, which component 
is U(l)pQ-noninvariant is just a matter of the choice: UIjq or U3jq is the lightest spectrum 
(u-quark). Here, we discuss only the case where UIjq is the lightest one. The resultant 
U(1)pq concerning the above constraints is 3 



Q 



10 



1. 



Ql = Ql 



h Q\ = o, 

Ql = Q 2 1= 2, Q\ = 2-q, 

Qm = -l - q, Qih = ~ 2 > 

where the superscript means the generation number. However, this assignment is not 
successful. For example, the lowest order term of the 11-component of lepton Yukawa 
coupling is (s^C/lsC/liL^g, and thus Qis®) = — <z/2- Determinant of lepton mass matrix is 
proportional to (s^) 2 , and so (s®) is 10 -3 M s t to fit the phenomenological values. Therefore, 
PQ symmetry breaking scale is around 10 15 GeV in this assignment. 

Up to now, we discussed the difficulty of obtaining the intermediate scale (~ 10 11 GeV) 
for the PQ symmetry breaking scale. We must resort to the GUT scale PQ symmetry 
breaking scale. For cosmological application, we may resort to the anthropic principle, 
a la Ref. p|. Then, not worrying about the PQ symmetry breaking scale, we look for 
approximate global symmetries surviving up to a high order. 



fields 


Qg\A 


fields 


QglA 


fields Q g iA 


fields QgiA 


fields 


QglA 


U3 5 


4 


Ul 5 


4 


T2 5 


4 


T4 5 (3) 


-4 


T6g(2) 


-4 


T6 5 (2) 


4 


C/3i 


-4 


Ult 


-4 


T2i 


-4 


T6i/(2) 


-4 


T6u(2) 


4 


^2oio 


2 


T2 + — 

Ol6 


1 


s u 


4 


-2(2) 


4 


-8(2) 


-4 


S ?(2) 


-4 


dr(2) 


4 


^(2) 


-4 


^( 2 ) 


-4 


C? 


4 


c° 


4 


c 7 ° 


-4 


ct 


4 


C+ 


4 


T1 1A 


-4 


Tli D 


4 


T1 1E 


4 


T1 XF 


4 


The 


4 


6i 


-4 


T1 + 


4 




-4 




-4 


T7 1A 


-4 


T7 1D 


4 


T7 1E 


4 


T7 1F 


4 


T7 1G 


4 


S3 


-4 


T7f D 


4 




-4 


T1- 1B 


-4 











Table 6: Non- vanishing U(l) g iA charges. Fields $('P) denotes the multiplicity T 5 in case V > 1. 

Up to .D = 7 superpotential terms, we have a global symmetry which we call U(l) g iA- 
This symmetry is flavor- independent, and so it belongs to the class described in Eqs. 
d4T3| ), ( gig ) and ( gig ). As explained, it does not have U(l) g i A -SU(5)-SU(5) anomaly. 
Low energy fields and their U(l) g iA global charges are listed in Table |(| From Table || and 



s Here, we assign Qj^ = 1 without loss of generality and Q\ = using U(i)x gauge symmetry. 
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|6|, the trace of the nonabelian anomalies are given by 



Tr (Qan?Su(5) r SU(5)) = ~9, 

" 9 ' (4.19) 



Tr (Q an Tso(lO)' r SO(10)', 



Tr (<3 g iA?SU(5) r SU(5)) - °> 
Tr (QgiA r SO (io)' r SO (io)') = 4 

We find that there is a gauge symmetry U(l) g iv which is very close to U(l) g iA, 

11111 

VglV = Z\ + 77-^2 + 77-^3 — + + 7f^6- 

OOOOO 



The U(l) g iv charges are given in ( plsD except one field T2j T5 for which Q glA (T2j T5 ) = -1. 

At the D = 8 superpotential terms, U(l) g iA is broken. The following superpotential 
terms at dimension 8 break U(l) g iA, 

AW = F(other fields)T2+_r2+ Tg T2 ol0 , (4.20) 

where 



Fi = T6 l0 T2 5 T4 5 T4 5 Cz + T6 10 T6 w Ti 5 U2- 5 C^ + T4 5 T6 5 U2- 5 T6~ 5 C^ 
+ T2 5 T4 5 T4- 5 T6- 5 Cz + T& W U2- 5 T4- 5 T6- 5 C 2 ~ + T6 10 T 4gT4gT2 1 C 2 - 
+ T6 10 U2- 5 TA- 5 C^T6 U + T6 W T6 W C$C%C 3 - + r2 5 T6gC 2 -sg 8 g 
+ T2 5 T6 b C 2 ~ s a 7 sl + T2 5 TQ l C^s^s^ + T2 b TQ l C^s^s^ 
+ T4 5 T45T2 1 C 2 -T6ij + TA^C^TGuTGu + TA h U2- 5 C^h 2 h x 
+ T4 5 U2 E C^h!h 2 + T4 5 U2 E C^h 4 h 3 + r4 5 C/25C 2 -/i3/i4 
+ r6 5 r6gC7 5 C 6 C73~ + T2 5 T6- 5 C° e h 3 Cz + T2 5 T6~ 5 Cfh 1 C^ 
+ T2 1 C~s 5 s° 6 T6 1 j + T2 1 C~s 7 s 8 T6 1J + T2 x C^s\s~ T6 U 
+ T2 1 C 2 ~s+SzT6 1 j + ClC^T6 u T6 u C^ + T2 1 C^T6 1J h 3 C^ 
+ T2 1 ClT6 u h 1 C 3 - (4.21) 



+ C*2 s ed 2 s 2 d 3 + C 2 s^s^d^d^ + C 2 Sg(i+ d% + C 2 s^ci^d^s^ 

+ C*2 S 5 S 6<^2~^2 + ^2 S 7 S 8^2~ ^2 + ^2 s 5 s 7^3 f ^2 + ^2 S 2~ ^2 S 2 
+ C*2 S 6 S 8<^2~^3 + ^2 S 5 S 6^3~^3 + ^2 s 7 a 8^3~^3 + ^2 s 2~ S 2 ^3 
+ C*2 ^2" S 3~ S 2 ^3 + ^2 S 5^2 S 2 ^3 + ^2 S 8 S 2 ^3 ^3 + ^2 s 2~^3~^2 S 3 

+ C 2 d 2 Sg"^ s 3 + C 2 ci 2 s 3 + C 2 d 3 s^dg s 3 + C 2 Sgd 2 o? 3 s 3 

(4.22) 
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(4.23) 



F 3 = C 6 °C 2 -sgs^2 + C%C 2 s 7 s° 8 h 2 + ClD 2 s+d 2 h 2 + C%C 2 s+s 2 h 2 

+ C%C 2 sp 3 h 2 + + C^D^s+d 2 hi + CgOf sjaj/u 

+ C^Df /i 4 + CgCjsga?^ + C$D 2 s+d 2 h 2 + ClD 2 s+d^h 2 
+ CgC^s+Sg + CgC^sgag^ + C 5 °C 2 -s?s^4 + ClC 2 s+s 2 h A 
+ ClD^s+d^hi + C^C 2 sjs^h 4 + C^^ 2 /iiC 3 _ + C^d+d 2 h 3 C^ 
+ C^d+d 3 h 3 C 3 - + C$C%h 2 h 3 C2 + C 5 fi+d 2 hiCJ + Cgd^fciGf 
+ ClClh 2 h x Cz + Cfcg/iifaGf + C 6 °4d3 /j 3 c 3 - + ClclhJi 3 Cz 

4.4 SM Yukawa couplings 

4.4.1 Q cm = I quark masses 

In the 10j5j5# couplings of Table ||, we can read off that the following singlet field must 
develop near GUT scale VEVs to have a realistic Q em = | mass matrix, 

huhuhs^ClCl (4.24) 

To achieve the hierarchy ^ ~ 1CP 6 among up-type quarks, some of the singlet vev might 
also have hierarchically small values ~ 10 -3 M st . 

4.4.2 Q cm = — i quark masses 

According to Table ||, the following singlets are required to develop near GUT scale VEVs 
toward a successful Q em = — \ mass matrix, 



hi,hi,h 3 ,h 3 . 



(4.25) 



These fields can have GUT scale VEVs and by some fine tuning we can obtain realistic 
Qem = — \ quark mass matrix. 

4.4.3 Q cm = — 1 lepton masses 

Finally consider the charged lepton mass matrix. According 5 i lj&H couplings of Table 
the following fields are likely to have some VEVs of order Mgut ~ Mj: 



hi,hx,h 3 ,h 3 , 



(4.26) 



Out of these, s 2 is required not to obtain a GUT scale VEV but a somewhat smaller 
scale at Mli- The GUT scale VEVs themselves may have a small hierarchy to fit to the 
experimental values for masses and mixing angles. The reason for choosing a smaller VEV 
for s 2 is the following. The charged lepton matrix 5jlj5# from Table |3] is proportional to 
(5h) times 



/o 

h 3 h 3 1 





hih 3 + hji 3 
h\h\ 



+ 



s 2 s 2 s 2 h 3 s 2 
s° 2 h 3 
y s$hi 



(4.27) 
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where we separate the mass matrix into the leading term and the perturbation in powers 
of s 2 - Eigenvectors of the left s 2 independent matrix are (1,0, 0) T with eigenvalue zero 4 
and the others which has the component only in the second and the third elements. This 
implies that the change of eigenvalues for the second and the third eigenstates are vanishing 
in the first order perturbation. Three eigenvalues in this order are given by 

m e = {slf 

m^ = \ (h 3 h 3 + hih - ^{h 3 h 3 - h^f + 4(1 + hJT 3 + /ii/13) 2 ) (4.28) 



m T = I [h 3 h 3 + h\hi + \J (h 3 h 3 - hihi) 2 + 4(1 + hih 3 + h x h 3 y 



However, for m e , the second order perturbation will also give (s 2 ) 2 order contribution, so 
we have to consider 0((s 2 ) 2 ) for an accurate estimate of electron mass. 

4.4.4 fi term 

Finally, let us examine the //-term among the relevant low energy couplings. Since the 
H parameter must be of order the electroweak scale which is negligible compared to the 
string scale, the \x parameter is presumably protected by a global symmetry J35[]. Under 
U(l) an and U(l) g iA, T\^U2^ is not invariant, and hence the ^-term is generated by some 
symmetry breaking VEVs of U(l) an and U(l) g iA- Up to dimension 7, we find that the /i 
term related couplings are the following 

' s° 2 sl, s u s 2 s° 2 , s^sKhsh + h 1 h 3 ),s u s 2 s° 2 (h 3 h 1 + hji 3 ), ' 
D+D^sl{hJi A + hJii),D+D 2 sl(h 2 h 3 + h 3 h 2 ), 
T4 5 U2s x { C$s%(sls%hx + d^d^hx + s^hi + d+d^hx), \ (4.29) 
Cis 2 s 2 slh 3 , Cls lS ° 4 s 8 h 3 , C§44<%h3, 
Cls^s^Cls^sfh^Cls^d+d^h 

These fi terms satisfy the U(l) an symmetry. Since TA^U2^{\.e.H u H ( i) carries 8 units of 
the U(l) an charge, the multiplied factor carries -8 units of U(l) an charge. If the singlet 
fields multiplied are flat directions, we can determine their VEVs by the minimization. For 
simplicity, we discuss the first two terms. Since the VEV of TA^U2i is nonzero, the factor 
s 2 s\ + s u (s 2 ) 2 must choose the minimum. If the VEVs of s 2 , s® and s u is not determined, 
probably being flat directions, this fi term determines a relation between s 2 ,s® and s u 
through the minimization condition. It is s 2 = —s\/2s u with the result fi = — (s®) 2 / 4s u . 
Since the electron mass is ~ (s 2 ) 2 , we take s 2 ~ 10 15 GeV. Then, the electroweak scale fi 
is obtained for example for s® ~ 10 5 GeV with an intermediate scale s u ~ 10 12 GeV. Of 
course, if they are not flat directions, their VEVs are determined by other more important 
terms. 

4.5 Search for an approximate symmetry 

For the sake of showing the existencee of vacua, let us choose only two kinds of VEVs, one 
kind carrying only the U(l) an and the other kind carrying only the U(l) g iA charge. In this 
way, we can separate two axion scales. 

4 The state (1, 0, 0) T will be interpreted as the electron. 
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From the Yukawa couplings and mass terms for the vectorlike pairs, we require near- 
GUT-scale VEVs to the fields: 

fci(0,0), hi (o,o), h 3 (0,0), h{o,o), 
C 3 °(-3,0), Cg(4,0), C 6 °(4,0), 
sg(-4 } 0), sr(-6,0), 
(4(5,0) or 4(5,0)) , 
(s?(-4,0) or S °(-4,0)), (4.30) 

where (Q an ,Q g iA) charges are shown. VEVs of these fields break U(l) an only. To break 
U(l)giA) we can consider the following 

S "(0,4), S 2(-4,4), S g(l,-4). (4.31) 

Thus, the domain wall numbers of the U(l) an and U(l) g iA axions are 9 and 1, respectively. 



5. Axion-Photon-Photon Coupling 

The PQ mechanism employs a global symmetry. The R symmetry is not considered for a 
Peccei-Quinn mechanism. This is because if we try to embed the axion field in the phase of 
a scalar field then this phase does not carry a SUSY parameter. So for the PQ mechanism 
gauginos are neutral. 

Let us consider the QCD axion and one hidden sector axion corresponding to SO(10)'. 
When we have two axions, we must consider the axion mixing and the cross theorem 
|2S, 14]. The higher axion potential corresponds to a lower axion decay constant. 



As discussed in the previous section, we have GUT scale VEVs for the fields of ( 4.3C| ) 



and intermediate scale VEVs Mj for the fields of ( 4.31| ). In fact, s° an d some other fields 



which appear in up-type Yukawa couplings have the VEV scale Mli somewhat larger than 
Mj but smaller than Mqut f° r a successful electron mass and for a ratio of up and top 
quark masses. The intermediate scale VEVs are in the region ~ 10 11 GeV. So, we study 
the following case, 

F A ~ 10 16 GeV, F G ~ 10 n GeV (5.1) 

where Fa and Fq correspond to axion decay constants derived from spontaneous symmetry 
breaking of U(l) an and U(l) g iA, respectively. It turns out that the decay constant of QCD 
axion is made to be of order GUT scale. 

The confining group we are interested in is SU(3) c x SO(10)'. We assume that a sub- 
group of SO(10)' such as SU(4) or SO(8) confines around lO 11 " 13 GeV. 

The early example of the axion mixing is between the MI axion and the MD breathing 
mode axion coupling as p9f. 



' -a A (FF + F'F') + oo 1 — a G (FF-F'F'). 



32vr 2 F M i v ' 327r 2 F MD 

But in our case the couplings are more complicated. The phases of C§(-3, 0), Cg(4, 0), 
Cg (4, 0), and s u (0, 4), s°(— 4, 4), Sg(l, —4), for example couple according to their U(l) an and 
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U(l) g iA charges. The axion decay constant can be calculated following the discussion of 
Sec. ^. Most probably, both axions, the QCD axion and the hidden sector axion, have the 
decay constants at the GUT scale because of the constraints on Yukawa couplings. 

Here, for the simplicity of illustration with relatively complex U(l) an and U(l) g iA 
charges, let us choose VEVs of (Cg(4, 0)) = V\ and (s°(— 4, 4)} = V 2 for global symmetry 
breaking. The relation s\ ~ — 2s u s\ of Subsubsec. 4.4.4] has been an illustrative example 



for flat directions, and here we present the coupling calculation to show the validity of 
cross theorem and the calculational method of axion-photon-photon coupling with a field 
having both global charges. So the VEV scale of s\ is chosen at an arbitrary scale here. 
The U(l) an breaking direction with Cg and s\ is 



l&V? + 4 2 V 2 2 = V (5.2) 
and the U(l) g iA breaking direction is 



I/PV* = AV 2 . (5.3) 

The normalized phase component for the U(l) an symmetry is 

W x a c - W 2 a s . 

O an = — ; = COS Oa-nCLn — Sin^anB.. 5.4) 

y/4?V? + 4 2 l/ 2 2 
The component orthogonal to a an can be written as 

a rtho = sin 6> an a c + cos 6> an a s , (5.5) 

or ac and a s are given by 

a c = cos # an a an + sin 6 an a OTtho , a s = - sin 6> an a an + cos 0an«ortho- (5-6) 
In the V± 3> V 2 limit, we have 

V 2 

cosfl an ~ 1, sin6U ~ — (5.7) 

Vi 

The phase component for the U(l) g iA symmetry is 

QgiA = the phase of s® = a s . (5-8) 

Since U(l) an -SU(5)-SU(5) and U(l) an -SO(10)'-SO(10)' anomalies are -9 units, and 
U(l)giA-SU(5)-SU(5) is and U(l) glA ~SO(10)'-SO(10) / anomaly is 4, we obtain the fol- 
lowing couplings 

.9^ -^(FF + F'F') + ( — ^r(F'F') 
V J 32tt 2 V J \ AV 2 J 32tt 2 V ' 



_L (-9^.) FF + — (-9^ + F'F' 
32vr 2 V V J 32vr 2 V V W 2 ) ' 



(5.9) 
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from which we obtain a potential proportional to 

- ~ a qcd cos f^ 5 c) - ^h m cos (jr 5 ^ ~ ( 5 - 10 ) 



where we used 



and took the limit 



*l = y, F 2 = ^j, ac = 4a c , a s = 4a s (5.11) 



Thus, the hidden sector axion decay constant is 



Vo 

F h = -f. (5.12) 



/.From Eq. the axion mass matrix is estimated as 



M^I^I^'^f] (5-13) 

where we simplified the confining scales as Aqcd and A^, and m is the mass scale of the 
hidden sector quarks. Introducing hierarchically small parameters 

e = ^2' r l = — la - * ( 5 - 14 ) 
we obtain the larger hidden sector axion mass and the smaller QCD axion mass 

2 A 4 (5-15) 

2 ^ A QCP 

m a,QCD — _p2 

which shows the validity of the cross theorem on the axion masses and decay constants. 

Until now in this section, we treat J7(l) g iA as an exact global symmetry. As mentioned, 
this symmetry is broken at order dimension 8 in the superpotential and thus we have to 
discuss the effect of its breaking on the axion. With the symmetry breaking terms, the 
vacuum may be shifted to a CP-violating one. The potential becomes 

.4 /a an \ 3 / a an a g i A \ {5m) 2 2 

- a qcd cos I — 1 - A^mcos I — + — 1 + — ^— (a g i A - flbr^i) , (5.16) 

where we denote the magnitude of the symmetry breaking term by (5m) 2 and the shift 
of the vacuum angle by #br- In the approximation that the shift is small, we linearize 
dV/da = 0. Then, 



A 4 , mkl w \ / 




Ajm Aim , , r \9. I I „ , . J I f,T™\2 
FiF 2 ~Ff" 



(5m)^ brJ F 2 



(5.17) 
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We have two cases : (i) 5m » -52-, (ii) dm <C -fe-- For the case (i), the QCD vacuum 
angle shift is 

^ ~ -flta. (5.18) 

So this cannot be used as a QCD axion. But if it belongs to the case (ii), we have 

a an (^m) 2 



p 1 A 3 r 



Ar (5.19) 



which can be sufficiently small. 

Let us now proceed to calculate the axion-photon-photon coupling. The global anomaly 

is 

U(l) an -^m-^T: -15 (5.20) 
U(l) glA -iC-^T: (5.21) 

We choose U(l) an and U(l) g iA for the axions. For illustration we present again for the case 
of C% and s% 

For the electromagnetic field Ftm, we have the following anomaly 
-15a an - 

^em^em 



32vr 2 V V 



1 ( -15(4a c + f y/ea s ) \ „ ~ 1 f -15a c 



32vr 2 I 4Vi / " n Lm ~ 32vr 2 V Vy 



^em-fern — on o I Jr I fem^em (5.22) 



1 



5 oqcd\ 



^em^em (5.23) 



32tt 2 \3 Fi ; 

whence we obtain the following high energy value for the axion-photon-photon coupling, 

qqcp / e 2 ~ \ 

^a,QCD V-327T Z / 

where 

Ca 77 = 3, F a ,QGD = Fi. (5.25) 

The axion-photon-photon coupling at low energy takes into account the QCD chiral sym- 
metry breaking and we obtain 



C(277 C<J77 



c xSB - -0.26 (5.26) 



where we take c x sb — 1-93. The value c x sb depends on the current quark mass ratio 
Z = m u /md, the instanton contribution factor £ = m s /A instanton , and r^ s = m^/m s . For 
£ = 0.1 and rds = ^j, we obtain 

|(4 + 1.05Z) 

C ^ = 1 + 1.05/ < 5 ' 27 > 
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Figure 2: Experimental bound from various experiments, especially including CAST 2003 data. 
This figure is taken from the result paper of CAST experiment Here, we show the prediction 
of Z12 model by the thick line. The blue band is the 20 % theoretical error of Rcf. 44|. 



Inclusion of one-loop effects and second order chiral perturbation give sizable contributions 
to squared meson masses |pf |. For Z = |, we obtain c x sb = 1-93. With the next order 
effect up to the 20% level, we obtain a band for c x sb, i-e. c x sb = 1-76 ~ 2.26. Thus, 
c a77 turns out to be sizable, which can be detected. Fig. ^ shows the current experimental 
search limit on c aTy and the prediction point of the present model with the band shown. 
The solid line corresponds to c x sb = 1.93. It is the first reliable calculation of the axion- 
photon-photon coupling. The cavity detectors [45] and the high Z Rydberg atom Kyoto 
axion detector Jig] already used this kind of axion-photon-photon interaction. But the 
recent CAST detector [^] seems to be the most promising one for detection of a very light 
axion in the region F a ~ 10 11 GeV. The detectability of axion using the c a77 coupling was 
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proposed by Sikivie [47]. However, the magnitude of F ,QCD in the model we studied is at 
the GUT scale and it cannot be detected in these kinds of detectors. 

What will be the case if we use VEVs of all the fields instead of just C® and S4? 
We have discussed at length that it was not possible to find a global symmetry which 
has the QCD anomaly consistently with the breaking scheme of the flipped SU(5). More 
importantly, the hidden sector axion potential is higher than the QCD axion potential and 
the cross theorem dictates that the decay constant of the QCD axion is the breaking scale 
of U(l) an . Thus, in our model the QCD axion cannot be made detectable. 

Finally, we mention that the hidden sector axion potential can be made smaller than 
the QCD axion potential by making the hidden sector quark mass extremely small in which 
case the cross theorem acts in the other direction. In this case, we can expect that the 
QCD axion can fall in the detectable region. However, we have not found any approximate 
global symmetry possessing SU(5) anomaly and hence this desirable scenario is not realized 
in the present model. 



6. Conclusion 

In this paper, we presented a general method to house a QCD axion in string-derived 
MSSM models. One related objective is to make it observable in ongoing or future axion 
search experiments since the axion derived from superstring might be the most significant 
prediction of string. We presented the criteria that should be satisfied in these models. 
Since the magnitude of the axion decay constant is essential in the solution of the strong 
CP problem, cosmology, and in axion search experiments, we presented a general formula 
for the decay constant in Eqs. (|2.7| ) and (2.9). We must consider the full Yukawa coupling 



structure of matter fields toward this objective. So far, this kind of full Yukawa coupling 
structure has not been studied except in a recent Z12-/ model |26[] . Here the Yukawa 
coupling has been given completely, which made it possible for us to pick up an approximate 
global symmetry so that a phenomenologically allowed QCD axion results. However, in the 
Z12-7 model we study there does not exists a vacuum where F a ~ 10 11 GeV with successful 
MSSM phenomenology. This might be a most probable situation in string models. It will 
be very interesting if one can find a string model with an observable QCD axion. 
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